Search for: All records

A preservative scheme is presented and analyzed for the solution of a quenching type convective-diffusion problem modeled through one-sided Riemann-Liouville space-fractional derivatives. Properly weighted Grünwald formulas are employed for the discretization of the fractional derivative. A forward difference approximation is considered in the approximation of the convective term of the nonlinear equation. Temporal steps are optimized via an asymptotic arc-length monitoring mechanism till the quenching point. Under suitable constraints on spatial-temporal discretization steps, the monotonicity, positivity preservations of the numerical solution and numerical stability of the scheme are proved. Three numerical experiments are designed to demonstrate and simulate key characteristics of the semi-adaptive scheme constructed, including critical length, quenching time and quenching location of the fractional quenching phenomena formulated through the one-sided space-fractional convective-diffusion initial-boundary value problem. Effects of the convective function to quenching are discussed. Numerical estimates of the order of convergence are obtained. Computational results obtained are carefully compared with those acquired from conventional integer order quenching convection-diffusion problems for validating anticipated accuracy. The experiments have demonstrated expected accuracy and feasibility of the new method. 

The aims of this paper are to investigate and propose a numerical approximation for a quenching type diffusion problem associated with a two-sided Riemann-Liouville space- fractional derivative. The approach adopts weighted Grünwald formulas for suitable spatial discretization. An implicit Crank-Nicolson scheme combined with adaptive technology is then implemented for a temporal integration. Monotonicity, positivity preservation and linearized stability are proved under suitable constraints on spatial and temporal discretization parameters. Two specially designed simulation experiments are presented for illustrating and outreaching properties of the numerical method constructed. Connections between the two-sided fractional differential operator and critical values including quenching time, critical length and quenching location are investigated.